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Abstract 

Four highly ordered hydrogen-bonded models of P- cyclodextrin (P-CD) and its inclusion complex with benzene were investigated 
by three different theoretical methods: classical quantum mechanics (QM) on AMI and on the BP/TZVP-DISP3 level of approxi- 
mation, and thirdly by classical molecular dynamics simulations (MD) at different temperatures (120 K and 273 to 300 K). The 
hydrogen bonds at the larger 02/03 rim of empty P-CDs prefer the right-hand orientation, e.g., 03-H 02-H in the same glucose 
unit and bifurcated towards • 04 and 03 of the next glucose unit on the right side. On AMI level the complex energy was 
-2.75 kcal moP 1 when the benzene molecule was located parallel inside the P-CD cavity and -2.46 kcal moP 1 when it was posi- 
tioned vertically. The AMI HOMO/LUMO gap of the empty P-CD with about 12 eV is lowered to about 10 eV in the complex, in 
agreement with data from the literature. AMI IR spectra displayed a splitting of the O-H frequencies of cyclodextrin upon com- 
plex formation. At the BP/TZVP-DISP3 level the parallel and vertical positions from the starting structures converged to a struc- 
ture where benzene assumes a more oblique position (-20.16 kcal moP 1 and -20.22 kcal moP 1 , resp.) as was reported in the litera- 
ture. The character of the COSMO-RS o-surface of P-CD was much more hydrophobic on its 06 rim than on its 02/03 side when 
all hydrogen bonds were arranged in a concerted mode. 

This static QM picture of the P-CD/benzene complex at 0 K was extended by MD simulations. At 120 K benzene was mobile but 
always stayed inside the cavity of P-CD. The trajectories at 273, 280, 290 and 300 K certainly no longer displayed the highly 
ordered hydrogen bonds of P-CD and benzene occupied many different positions inside the cavity, before it left the P-CD finally at 
its 02/03 side. 
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Introduction 

Cyclodextrins (CD) are a family of conical shaped cyclic 
oligosacharides consisting of 6-8 (or up to 10) glucopyranose 
units linked by a-(l— >4) glycosidic bonds. At the narrower rim 
of the truncated cone (06 side) there is one primary hydroxy 
group per glucose unit whereas at the wider rim there are two 
secondary hydroxy groups (02/03 side) [1]. Cyclodextrins have 
many possibilities of forming hydrogen bonds [2]: firstly, 
cyclodextrin monomers form some intramolecular hydrogen 
bonds or closed rings of hydrogen bonds at both rims of the 
cone; secondly, they form intermolecular hydrogen bonds with 
water molecules in aqueous solution or with included guest 
molecules that fit into their cavity. This fit is extremely specific: 
the orientation of one and the same guest molecule inside the 
cyclodextrin cavity may be opposite in solution and in crys- 
talline state, as was found for 4-fluorophenol in a-CD [3]. The 
reactivity of aromatic guest molecules, radicals or excited 
states, was found to be altered because of complex formation 
with cyclodextrins [4]. 

Cyclodextrins form three types of dimers, 02/03 to 02/03, 02/ 
03 to 06, and 06 to 06. They can also associate to extended 
stacks in crystalline state or solution, with guest and solvent 



molecules located inside the cavities, between dimers, or in the 
channels next to the stacks [5]. The influence of cyclodextrin on 
a broad variety of molecules was reported: from changes of a-/ 
P-sheet ratio upon complex formation of cyclodextrins with side 
chains of proteins [6], to size control of the electrostatic self 
assembly of nanoparticles [7]. Also a pH-controlled release of 
many cyclodextrins in long stacks on polymers, such as polyro- 
taxanes, combined with mesoporous silica particles was 
observed [8]. 

Complex formation with solvent molecules can take place at the 
same time in different positions: the crystal structure of 
heptakis(6-0-triisopropylsilyl)-P-CD benzene pyrene solvate is 
a nice example for competing solvents that intercalate in the 
cavity of one P-CD (benzene), between P-CD dimers (pyrene), 
and in channels outside P-CDs (benzene). The P-CD dimers are 
connected with hydrogen bonds on their 02/03 sides 
(Figure 1). 

Investigations of such "flat energy hypersurfaces" due to 
hydrogen bonds and noncovalent interactions demand several 
theoretical methods to capture the entire network of forces on 




Figure 1: Crystal structure of heptakis(6-0-triisopropylsilyl)-(3-cyclodextrin benzene pyrene solvate; [Cio5H2io035Si7,o.5(Ci6 H io)>3.5( c 6 H 6)L taken 
from The Cambridge Crystallographic Data Centre CCDC [9]. Colour code: benzene: red; pyrene: magenta; cyclodextrin corpus: light green; 
cyclodextrin side chains: dark green; hydrogen bonds: light blue, dotted lines. 
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which they rely. Molecular dynamics simulations displayed 
different hydrogen bond patterns of cyclodextrins in crystal and 
in solution and confirmed the experimental findings of soluble 
cyclodextrin complexes of cholesterol type versus precipitates 
of czs/£rafts-cyclohexadecenone//CDs [10-12]. MD with 
^-dynamics and calculation of the solvation-free-energy differ- 
ences was used to distinguish amylose helices from cellulose 
sheets by analysing the different reactivity of oxygen atoms 02, 
03 and 06 of the sugar units with and without methylation, in 
line with experimental data [13]. Cyclodextrin-complex forma- 
tion with substituted benzenes shows multiple topologies/ 
configurations (guest up/down) in MD with ^-dynamics [14]. 
The association constant K a for a- and P-CD inclusion 
complexes with several benzene derivatives was investigated by 
a genetic algorithm [15] and neuronal networks [16]. An 
overview of quantum mechanical methods to study cyclodex- 
trin chemistry is given in [17]. The COSMO-RS solvation 
model [18] allows for calculations of several thermodynamic 
properties [19] once the polarisation charge surface of mole- 
cules has been determined by DFT calculations. With the 
current implementation of COSMO-RS, hydrogen-bond ener- 
gies are predicted within 0.36 kcal mol -1 relative to QM dimer 
calculations [20]. 

The necessity to understand the function of cyclodextrins in all 
details arises from their many applications in industry [21], 
agriculture, food [22] and health care [23]. Cyclodextrins and 
their complexes are produced in large industrial processes [24- 
26]. Our intention for this study was to find a "modular, robust 
molecular-modelling workflow" for how to analyse the forma- 
tion of cyclodextrin complexes for many guest molecules in 
general, before extensive experimental work is carried out. 
Lemon grass oil [27] or retinol [23] can serve as typical molec- 
ular examples of the "practical industrial workbench". To reach 
reasonable structures of these complexes we need to construct 
several models because of (a) isomers, (b) n:m guest/host stoi- 
chiometry, (c) preferred cyclodextrin cavity size, etc. We start 
here with the first step, the empty P-CDs, their subtle but 
strongly influential hydrogen bonds, and their complex forma- 
tion with the simplest aromatic molecule of highest symmetry 
without substituents, i.e., benzene [28]. 

Methods 

While calculating several cyclodextrin complexes with molec- 
ular dynamics [29], semi-empirical AMI, classical quantum 
mechanics [30] and COSMO-RS (Turbomole) [31], we 
observed that the orientation of all hydrogen bonds of the 
cyclodextrins already played a crucial role when the initial 
structures were constructed. It was not sufficient just to take the 
"normal route" of using the energetically lowest-lying struc- 
tures from MD trajectories, and then proceed to AMI and 



further to QM optimizations; on the contrary, these structures 
turned out not to be the best ones possible. The best ones in QM 
finally were those that had hydrogen bonds oriented as symmet- 
rically as possible. Therefore, we started with such manually 
constructed models and their AMI optimisations. 

Step 1: Our four P-CD models are conformers; all the hydrogen 
bonds of each rim are oriented in the same direction. We named 
them BCD0231061, BCD023r061, BCD023106r and 
BCD023r06r, referring to their manually constructed 
hydrogen-bond orientations when looking either at their 02/03 
or 06 rim (1 = left hand, r = right hand orientation) by using 
Gauss View [30]. Each model was optimised until convergence 
was reached "AMI (opt=calcall, verytight)" by using the 
Gaussian9W program [30]. Molecular properties, such as 
energy, dipole moment, HOMO/LUMO molecular orbitals and 
IR spectral data, were analysed with GaussView. 

Step 2: The four AMI optimised P-CD models were used as 
starting structures in our QM calculations by using Turbomole 
v. 6.4 with COSMO-RS [31]. Two models were calculated for 
each structure, i.e., the structure in vacuo (.energy files) and the 
COSMO-RS structure (.cosmo files), the latter being optimised 
in a dielectric field with the dielectric constant of water (the 
value of the dielectric constant is a variable in the COSMO-RS 
program and can be chosen by the user). Up to now it was 
possible to calculate molecular structures with COSMO-RS on 
DMOL, BP/SVP or BP/TZVP level of approximation. Our trial 
to optimise the P-CD complexes with BP/TZVP failed: not 
always but in several cases, here, for example, the benzene 
guest was expelled for the BCDbenzene023106r conformers. 
Now, however, the current version of Turbomole allows the 
calculation of more detailed energy hypersurfaces especially for 
hydrogen bonds (keyword BP-TZVP-DISP3) since Grimme's 
dispersion corrections for nonbonded interactions were imple- 
mented recently [32]. The current release of COSMO-RS 
C30 1201 allows, for the first time, a more detailed o-surface to 
be calculated (BP/TZVPD-FINE), but for the time being as 
single points (SP) only. COSMO-RS structures are saved in 
COSMO databanks to perform graphical and thermodynamic 
analysis with COSMOthermX [19]. In this study we used the 
BP/TZVP-DISP3 method to calculate the P-CD models, and 
with this method no benzene molecule was expelled from the 
P-CD cavity any more. 

Step 3: Starting from the neutron diffraction study of P-CD [33] 
we constructed one model named "BCDcry" of one P-CD by 
adding hydrogen atoms using the Visualizer of Materials Studio 
v. 5.5 [29]. Additionally, we used the four empty P-CDs of the 
in vacuo runs (BP/TZVP optimised .energy files) and the two 
best structures from the AMI optimisations, BCD023r061/ 
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benzene parallel and vertical. With Materials Studio molecular 
dynamics (Forcite plus program and Compass force field) we 
simulated their trajectories at 300 K. In detail: After force-field 
geometry optimisations the structures were heated up each from 
0 to 3 K, to 100 K, to 200 K, and to 300 K. Finally at 300 K a 
trajectory of 6000 ps length was calculated using an integration 
step of 1 fs. From the last 2000 ps of the trajectories the aver- 
aged total energy E tot <4000-6000>, averaged hydrogen-bond 
energy £ HB <4000-6000>, and the averaged hydrogen-bond 
distance and angle were calculated. 

Next, the energetically two best optimised AMI structures 
(BCD023r061 with benzene included in parallel and vertical 
positions) were used as starting structures and simulated at 300, 
290, 280, and 273 K for 500 ps each, and their release of the 
guest was analysed from these trajectories. 

Results and Discussion 

AM1 Calculations 

The four conformers of empty p-CD 

The BCD023r061 structure turned out to be the lowest energy 
conformer, also having the lowest dipole moment of only 
0.5618 Debye (dipole-moment components x = 0.0; y = 0.0; 
—z = 0.5618) pointing from the 06 rim towards the 023 rim of 
the cone. This structure is so symmetric that all its Mulliken 
partial atomic charges on all types of oxygen atoms were iden- 
tical, e.g., equal to each average number (Table 1). Second, with 
about 6 kcal mol -1 above, was conformer BCD023r06r. This 
oval structure has the highest dipole moment of 2.5639 D (x = 



+0.0268; y = +0.0236; z = -2.5637) and its 06 atomic charges 
display a range from -0.340 to -0.349 (average -0.344), but all 
02 atoms still have the same charge of 0.325, and 03 and 04 
oxygen charges vary very little (-0.337 to -0.338 and -0.300 to 
-0.301, respectively). The third conformer BCD0231061, lying 
about 7 kcal mol -1 above, was again total symmetric as shown 
by its dipole moment components of x = 0.0; y = 0.0; z = 
-1.0687 and also with respect to all its oxygen charge values. 
The fourth conformer, BCD023106r, had the highest energy of 
about 12.7 kcal mol -1 compared to the ground state 
BCD023r061 and the second largest dipole value of 2.4223 D 
(jc = +0.5806; y = -0.0872; z = -2.3501). Its 06 oxygen atoms 
display a broad range of Mulliken partial charges from -0.338 
to -0.360, although their average hardly differs from the 
average values of the other 06 types (see Table 1); 02 and 03 
have small ranges (-0.333 to -0.335 and -0.345 to -0.347) but 
some differences can be observed in the 04 values, i.e., -0.288 
to -0.296. These subtle differences were just given as a few 
examples here in the text (e.g., they are not given in Table 1) 
since they may be of importance during guest inclusion or re- 
activity. The Mulliken partial atomic charges of the benzene 
molecule in symmetry were C: -0.130 and H: +0.130. 

The four conformers of the p-CD complexes with 
benzene in a parallel or vertical position 

All four conformers of the complexes preferred the inclusion of 
benzene in a parallel position (Table 2). The energetic order of 
the benzene complexes and their sequence of dipole moments 
was qualitatively the same as for the empty ones: BCD023r061, 



Table 1: AM1 results for the four empty [3-CD conformers. AM1 optimised energies [a.u.], relative energies [kcal mol -1 ], dipole moments [D], LUMO/ 
HOMO energies [a.u.] and gap [eV], and averaged Mulliken partial atomic charges of oxygen atoms. Benzene molecule with D§ h symmetry in vacuo. 
AE is the energy of each conformer with respect to the lowest conformer in this table, which is 023r06l (bold). 



AM1 

p-CD 

empty 3 


abs. E 
[a.u.] 


AE 

[kcal mol -1 ] 


M 
[D] 


LUMO/HOMO 
[a.u.]; 
gap [eV] 


02 
atomic 
charge 


03 
atomic 
charge 


04 
atomic 
charge 


06 
atomic 
charge 


023I06I 


-2.653658 


7.06 


1.0687 


+0.06271 
-0.38469 
12.17 


-0.335 


-0.347 


-0.291 


-0.348 


023l06r 


-2.644708 


12.68 


2.4223 


+0.06083 
-0.38424 
12.11 


-0.334 


-0.346 


-0.293 


-0.346 


023r06l 


-2.664914 


0.00 


0.5618 


+0.06072 
-0.38650 
12.17 


-0.325 


-0.337 


-0.297 


-0.348 


023r06r 


-2.655373 


5.99 


2.5639 


+0.05835 
-0.38786 
12.14 


-0.325 


-0.337 


-0.301 


-0.344 


benzene b 


+0.034953 




0.0000 


+0.02039 
-0.35474 
10.21 











a AII four [3-CD conformers reached a RMS gradient norm of 0.00000001 [a.u.] in the geometry optimisations. 
b Benzene reached a RMS gradient norm of 0.00000002 [a.u.]. 
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Table 2: AM1 results for the four (3-CD conformers with benzene placed in a parallel or vertical position. AM1 absolute energies, relative energies and 
dipole moments for the complexes: columns 1-3; for benzene alone: columns 4 and 5. 

AM1 complex complex complex benzene benzene 
(3-CD/benzene complex 3 abs. E[a.u.] AE [kcal mol -1 ] |j [D] abs. E[a.u.] m [D] 


023I06I 












parallel 


-2.62313280 


7.0374744 


0.9544 


0.03501077 


0.0446 


vertical 


-2.62266055 


7.3338160 


0.7084 


0.03499769 


0.0092 


023l06r 












parallel 


-2.61493920 


12.1790404 


2.9750 


0.03499519 


0.0174 


vertical 


-2.61452474 


12.4391182 


1.5541 


0.03500111 


0.0104 


023r06l 












parallel 


-2.63434772 


0.00000000 


0.4399 


0.03502037 


0.0543 


vertical 


-2.63388503 


0.29034260 


0.1579 


0.03499812 


0.0096 


023r06r 












parallel 


-2.62574636 


5.39743941 


2.5365 


0.03499903 


0.0256 


vertical 


-2.62459809 


6.11799032 


1.9127 


0.03499288 


0.0138 


a AII four complexes reached a RMS gradient norm < 3 x 10 


3 [a.u.] in the geometry optimisations. 







BCD023r06r, BCD0231061, and BCD023106r. Benzene was 
slightly distorted in each cavity, as can be seen from its energy 
and dipole moment inside the complex. Inside the P-CD/ 
benzene complex the benzene neither in a parallel position nor 
in vertical position adopted an exact central position as 
displayed in Figure 2. 



Our calculated AMI HOMO/LUMO gaps of empty P-CD 
conformers, of benzene (Table 1) and the eight complexes 
(Table 3) are of the same order of magnitude as the values 
calculated with PM3 in the literature [34] where they found for 
P-CD 12.44 eV, for benzene 10.15 eV and for the P-CD/ 
benzene complex 10.15 eV. These authors found P-CD/benzene 




Figure 2: The two AM 1 -optimised stable conformers of BCD023r06l with benzene in a parallel (a and b) and vertical position (c and d). Top view 
from the 06 side. 
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Table 3: AM1 results for the four (3-CD/benzene complex conformers with benzene placed in a parallel or vertical position. LUMO/HOMO energies 
and band gaps: columns 1-3. Absolute energies of complex formation AE C ompiexation : columns 4,5. Data of column 5 are given in kilojoules per mole 
(kJ mor 1 ) for easy comparison to the text above and [34]. 

AM1 HOMO LUMO HOMO/LUMO gap 

P-CD/benzene [a.U.] [a.U.] [eV] A^complexation A^complexation 

complex [kcal mol 1 ] [kJ mol 1 ] 



023I06I 

parallel 

vertical 

023l06r 

parallel 

vertical 

023r06l 
parallel 
vertical 

023r06r 

parallel 

vertical 



-0.34409 
-0.34163 

-0.35373 
-0.34592 

-0.34784 
-0.34530 

-0.35776 
-0.35353 



0.03063 
0.03241 

0.02092 
0.02861 

0.02688 
0.02866 

0.01689 
0.02046 



10.196 
10.178 

10.195 
10.192 

10.197 
10.176 

10.195 
10.177 



4.29 
4.58 

9.43 
9.69 

-2.75 
-2.46 

2.65 
3.37 



17.93 
19.17 

39.44 
40.53 

-11.52 
-10.30 

11.07 
14.08 



of 1 : 1 inclusion type in aqueous solution by experimental UV 
absorption spectra with various concentrations of P-CD. From 
their UV plots they calculated their stability constants K s at 
three temperatures to be 1 16 L mol -1 at 291 K, 86 L mol -1 at 
298 K and 68 L mol -1 at 303 K. According to the equation 
AG = -RT In ^ s they obtained the corresponding three free 
energy values AG of -11.5 kJ mol -1 , -11.0 kJ mol -1 and 
-10.6 kJ mol -1 , indicating rather weak interactions between the 
guest and host molecule, but still a thermodynamically 
exothermic reaction of complex formation. 

Here, our best conformer, BCDbenzeneparallel023r061, 
displays an energy of complex formation 

"^complexation — ^complex ~~ ^BCDempty — ^benz/) 6 ^ (1) 



of -2.75 kcal mol -1 (Table 3), followed by the complex 
conformer BCDbenzenevertical023r061 with -2.46 kcal mol -1 . 
This is in very good agreement with the experimental value 
from [34]. All our other complex conformers have positive 
A^complexation values, up to BCDbenzenevertical023106r with 
+9.69 kcal mol -1 , and thus are endothermic, Table 3. 

IR spectra of empty p-CD conformers and their 
inclusion complexes 

The most important lines of the calculated AMI IR spectra are 
given in Table 4. The IR bands of benzene at 744, 1 145, 1579 
and 3194 cm -1 were characteristically shifted in the parallel 
versus the vertical complex. Especially the 3194 cm -1 H-C 
benzene stretch was split into a range from 3150 to 3199 cm -1 . 
Very interesting and easy to see was the concerted movement of 
all hydrogen bonds at the 023 and 06 rim: the empty P-CD's 



Table 4: AM1 results for benzene and the two stable 023r06l (3-CD conformers with benzene placed in a parallel or vertical position. All IR bands 
[cm -1 ] and relative intensities are listed for benzene; in the case of (3-CD, single bands are included only if their relative intensity is >200 or for "very 
close frequency clusters" (bold) 3 . 



AM1 rel. intensity rel. intensity rel. Intensity rel. Intensity comment 

IRfreq. benzene P-CD P-CD P-CD 

[cm -1 ] D 6/7 sym. 023r06l 023r06l 023r06l 

empty parallel vertical 

443 H(06) up/d 

444 729 713 H (06) up/d 

445 741 H(06) up/d 



511 276 H(02)up/d 

513 2 x 346 286 H (02) up/d 

514 240 H(02)up/d 
516 335 230 H (02) up/d 
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Table 4: AM1 results for benzene and the two stable 023r06l (3-CD conformers with benzene placed in a parallel or vertical position. All IR bands 
[cm -1 ] and relative intensities are listed for benzene; in the case of (3-CD, single bands are included only if their relative intensity is >200 or for "very 




744 


127 




up nn/rl 
n — \~> u|j/u 


757 




122 


H-C up/d 


760 


73 




H-C up/d 


1145 


2 x 1.08 




H-C in plane 


1 148 




I .0£- 


M P in nlanp 

1 1 w 1 1 1 \J\C\ \ \ c 


1150 


1.74 




H-C in plane 


1150 


1.51 




H-C in plane 


1152 




0.54 


H-C in plane 


1346 




344 


04 and 05 


1347 


471 491 




04 and 05 


1347 


471 355 




04 and 05 






i 


C~\A anrl 
vJf dllU vJO 


1351 


2 x 285 




C-06 frame 


1356-1357 


534(3) a 




C-frame 


1357 


2 x 312 




C-frame 


1357-1358 




562(3) a 


C-frame 


1418 


346 


338 


H-C up/d 


1419 


339 




H-C up/d 


1577 




14 


C=C 


1578 


15 




C=C 


1579 


2x12 15 


10 


C=C 


3150 




16 


H-C benz. 


3158 




22 


H-C benz. 


3165 




91 


H-C benz. 


3170 




35 


H-C benz. 


3171 


1D7 




H— C hpn7 
n — o uci iz.. 


3175 


2 x 88 75 




H-C benz. 


3186 


12 




H-C benz. 


3189 




20 


H-C benz. 


3194 






n — o uci iz.. 


3199 




17 


H-C benz. 


3414-3421 




879(7) a 


H-02 stretch 


3418-3421 


861 (7) a 




H-02 stretch 


3419 


2 x425 




H-02 stretch 


3447-3455 




733(5) a 


H-06 stretch 


3449-3451 


1004(7) a 




H-06 stretch 


3453 


2 x 536 




H-06 stretch 


3455-3456 


563(7) a 




H-03 stretch 


3457 


532(3) a 




H-03 stretch 


3456-3461 




832(9) a 


H-06/H-03 stretch 


a Sum of (a?) intensities within the given frequency range for a cluster of the same type (bold). 
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H-02 stretch at 3419 cm -1 was spread out to a range from 3414 
to 3421 cm -1 , the H-06 stretch at 3453 cm -1 was spread out 
from 3447 to 3453 cm" 1 , and the P-CD's H-03 stretch at 
3457 cm -1 from 3455 to 3456 cm -1 as a result of the interaction 
with the guest molecule. In general, the shifts were subtle; it 
seemed that the benzene itself was influenced very little if at all, 
the cyclodextrin on the contrary was substantially influenced in 
its frame vibrations over a wider range, its splitting most prob- 
ably being caused by the "rigid benzene wheel" it had included, 



see Figure 3. Interestingly, the concerted hydrogen-bond move- 
ments at the rims of P-CD rims were easily identified, for 
example in the BCDbenzeneparallel023r061 complex. Here, at 
445 cm -1 , all hydrogen atoms bound to 06 atoms concertedly 
moved up and down the almost perfect plane of 06 oxygen 
atoms (marked H(06) up/d in Table 4). 

Since our two stable P-CD/benzene complex conformers 
(parallel and vertical) were energetically close neighbours with 
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Figure 3: (a) AM1 calculated IR spectra of BCD023r06l empty and (b) the BCD023r06l/benzene inclusion complex with benzene in a parallel orien- 
tation. 
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only 0.29 kcal mol -1 energy difference calculated at 0 K, the 
experimental spectrum at room temperature should display all 
cyclodextrin absorptions, i.e., both our conformers exist at the 
same time. 

Summarised result of the AM1 calculations 

AMI energy/geometry optimisations in vacuo at 0 K showed 
that the (3-CD/benzene complex with all hydrogen bonds at the 
02/03 rim in a right-hand orientation and at the narrower 06 
rim in a left-hand orientation can include benzene in a parallel 
(best energy) or vertical position (higher in energy by only 
0.29 kcal mol -1 ); both conformers were thermodynamically 
allowed. The HOMO/LUMO gap of the empty P-CD with about 
12 eV was lowered to about 10 eV in the complex. The P-CD/ 
benzene inclusion complex was formed only by weak noncova- 
lent interactions, which influence the IR lines of the P-CD most 
at its frame frequencies and at highly ordered concerted move- 
ments of its hydrogen bonds. The IR frequencies of benzene 
were influenced only marginally upon inclusion. 

COSMO-RS calculations with Turbomole 
Quantum mechanical geometry/energy 
optimisations 

For the four empty P-CD models, for benzene in geometry, 
and for all eight BCD/benzene complexes, quantum mechan- 
ical geometry/energy optimisations with Turbomole and the 
implemented COSMO-RS method were performed at 0 K. For 
all starting structures the AMI optimised geometries were used, 
and the results are given in Table 5 for the calculations in 
vacuo, and in Table 6 for the calculations in the COSMO-RS 
dielectric field of water (no explicit water molecules, but 
instead dipoles that model the dielectric field of the solvent are 
used by this method [18]). 



The empty BCD023r061 model is still the lowest energy 
conformer among the four in vacuo (Table 5), but now 
the BCD0231061 is second and energetically close 
(0.26 kcal mol" 1 ), followed by BCD023r06r with 2.17 kcal 
mol" 1 and BCD023106r with 2.3 kcal mol" 1 higher. In aquo 
(Table 6) the empty hydrogen bond model BCD023r061 was 
again the most preferred (Figure 4). 

Although all optimisations were started from the AMI -opti- 
mised parallel and vertical geometries (from which the naming 
stems), the BP/TZVP-DISP3 method with its enhanced descrip- 
tion of dispersion forces led in general to structures in which the 
benzene molecule adopts a more oblique position [14,20,32]. 
The whole energy hypersurface of P-CD complexes is rather 
flat, and therefore it is difficult to identify the many local 
minima [17]. Therefore, geometry optimisations of P-CD 
complexes still remain a multiple-minimum problem in general. 
The COSMO-RS method combines these identified multiple- 
minima structures as a set of conformers in a databank (which is 
updated principally when new minima are identified by more 
sensitive theoretical methods). From these sets the thermody- 
namic properties of the material at various temperatures can be 
calculated (by analytical formulas with COSMOthermX) more 
realistically than for one structure only [18-20]. Molecules 
adopt many conformers at temperatures higher than 0 K, which 
will be shown for the P-CD/benzene complex by MD simula- 
tions in the next section. The energetically best structure of the 
BCD/benzene complex in vacuo due to the BP/TZVP-DISP3 
method resulted from the BCD023r061/benzenevertical starting 
structure. The AMI favourite BCD023r061/benzeneparallel 
model followed second with 0.06 kcal mol" 1 above, and the 
others were higher up to a maximum of 2.34 kcal mol" 1 
(BCD023106r/benzeneparallel), Table 5. All optimised BP/ 



Table 5: Optimised BP/TZVP-DISP3 energies in vacuo. Lowest energy conformer of each group in bold. 



BP/TZVP-DISP3 in vacuo 


£abs. [a.u.] 


E re |. [a.u.] 


E re | [kcal mol 1 ] 


BCD023I06I 


-4277.4641711793 


0.0004226810 


0.2652365729 


BCD023l06r 


-4277.4609298546 


0.0036640057 


2.2992002354 


BCD023r06l 


-4277.4645938603 


0.0000000000 


0.0000000000 


BCD023r06r 


-4277.4611238596 


0.0034700007 


2.1774601643 


Benzene 


-232.3358156576 


0.0071740456 


4.5017853670 


parallel023l06l 


-4509.8320598074 


0.0005758687 


0.3613633618 


parallel023l06r 


-4509.8289075499 


0.0037281263 


2.3394365033 


parallel023r06l 


-4509.8325394613 


0.0000962148 


0.0603757489 


parallel023r06r 


-4509.8291543743 


0.0034813018 


2.1845516926 


vertical023l06l 


-4509.8321043279 


0.0005313482 


0.3334263276 


vertical023l06r 


-4509.8289158257 


0.0037198504 


2.3342433180 


vertical023r06l 


-4509.8326356761 


0.0000000000 


0.0000000000 


vertical023r06r 


-4509.8291521426 


0.0034835335 


2.1859520815 
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Table 6: Optimised BP/TZVP-DISP3 COSMO-RS energies in aquo. Lowest energy conformer of each group in bold. 



or/ 1 ZVr-Ulbro in aquo 


£abs. [a-uj 


E re \. [a.u.] 


1 — ri . _ _ 1 — In 

t re i [kcal mol J 


BCD023I06I 


-4277.5397933200 


0.0044818408 


2.8123999202 


BCD023l06r 


-4277.5390925770 


0.0051825838 


3.2521231603 


BCD023r06l 


-4277.5442751608 


0.0000000000 


0.0000000000 


BCD023r06r 


-4277.5432219210 


0.0010532398 


0.6609185069 


Benzene 


-232.3408624099 


0.0021272933 


1.3348978187 


parallel023l06l 


-4509.9064609135 


0.0041674010 


2.6150858014 


parallel023l06r 


-4509.9056675976 


0.0049607169 


3.1128994621 


parallel023r06l 


-4509.9106283145 


0.0000000000 


0.0000000000 


parallel023r06r 


-4509.9096646704 


0.0009636441 


0.6046963092 


vertical023l06l 


-4509.9004752425 


0.0101530720 


6.3711542109 


vertical023l06r 


-4509.9055329919 


0.0050953226 


3.1973658847 


vertical023r06l 


-4509.9106053475 


0.0000229670 


0.0144120224 


vertical023r06r 


-4509.9095882112 


0.0010401033 


0.6526752221 




Figure 4: The lowest energy empty conformer BCD023r06l, COSMO-RS structure optimised with BP/TZVP-DISP3 (solvent = water), o-surface with 
charge code: oxygen red, hydrogen blue, hydrophobic area yellow and green. View a: from 02/03 side, view b: from 06 side, view c: from 02/03 
side. 



TZVP-DISP3 COSMO-RS energies in aquo were lower than 
the corresponding ones in vacuo, which is reasonable (Table 6). 
Here the BCD023r061 hydrogen-bond starting models 
remained the best ones after energy/geometry minimisations, 
for the empty P-CD and for the complex as well. 

The BP/TZVP-DISP3 energies of complex formation were 
calculated from the data of Table 5 and Table 6 by Equation 1 
and are shown in Table 7. All energies of complex formation 
A^Complexation in vacuo and in aquo were negative, indicating 
that the complexes will be formed. The complex structures of 
the reactions with the lowest energies were BCD023r061/ 
benzeneparallel in aquo and BCD023r061/benzenevertical in 
vacuo. 



Figure 4a-c and Figure 5a-c display different views of the best 
BCD023r061 empty model and the BCD023r061/benzene 
complex with their COSMO-RS o-surfaces. Hydrophilic areas 
are coloured in red (negative charge of oxygen atoms) or blue 
(positive charge of hydrogen atoms), which were mostly located 
at the 023 rim. Hydrophobic areas (coloured green/yellow) can 
be observed easily inside the BCD023r061 empty model 
(caused by its 06 highly ordered hydrogen-bond rim) and the 
whole surface of benzene while sitting inside the cyclodextrin 
cavity. 

Consequences for thermodynamic analysis 

The basis of thermodynamic analysis with COSMOthermX [18- 
20] is the so called o-profiles and o-potentials of each molecule, 
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Table 7: BP/TZVP-DISP3 energies of complex formation in vacuo and in aquo, calculated with the absolute data of BCD023r06l for EecDempty- Rela- 
tive AEcomplexation with respect to the energetically lowest complex conformer, column 2 and column 4. 



BP/TZVP-DISP3 


abs. AEcomplexation 
in vacuo 
[kcal mol" 1 ] 


rel. AEcomplexation 
in vacuo 
[kcal mol -1 ] 


abs. AEcomplexation 
in aquo 
[kcal mol -1 ] 


rel. AEcomplexation 
in aquo 
[kcal mol -1 ] 


parallel023l06l 


-19.861 


0.361 


-13.381 


2.615 


parallel023l06r 


-17.883 


2.339 


-12.883 


3.113 


parallel023r06l 


-20.162 


0.060 


-15.996 


0.000 


parallel023r06r 


-18.038 


2.184 


-15.391 


0.605 


vertical023l06l 


-19.889 


0.333 


-9.625 


6.371 


vertical023l06r 


-17.888 


2.334 


-12.798 


3.198 


vertical023r06l 


-20.222 


0.000 


-15.981 


0.015 


vertical023r06r 


-18.036 


2.186 


-15.343 


0.653 




Figure 5: The lowest energy complex conformer BCD023r06l/benzene; benzene occupies an oblique position inside the cavity. COSMO-RS struc- 
ture optimised with BP/TZVP-DISP3 (solvent = water), a: no surface, b,c: o-surface with charge code: oxygen red, hydrogen blue, hydrophobic area 
yellow and green. All views from 02/03 side. 



which show the characteristic differences in the electrostatic 
molecular surfaces that are calculated quantum mechanically at 
0 K for molecular conformers of up to 10 kcal mol -1 in the 
standard COSMO databases. These o-profiles and o-potentials 
were calculated with COSMOthermX for the empty (3-CD 
models and are displayed in Figure 6 and Figure 7. Especially 
the o-potentials (Figure 7) showed that the empty [3-CD models 
were split into two groups. One group had the two hydrogen- 
bond models with all hydrogen bonds in the left-hand orienta- 
tion (0231), the other one with the two models with right-hand 
(023r) orientation. This mirrored quantitatively the result of 
how much the two energetically preferred 023r models were 
closer to each other and were most different from the two 0231 
models, regardless of which orientation the hydrogen bonds had 
at the 06 rim. By the time these four empty [3-CD models have 
been combined into one set of [3-CD conformers and stored like 
this in the COSMO database, thermodynamic analysis with 
COSMOthermX will take into account these differences and 
give more accurate results for many material properties over a 
wider temperature range than with simple "one-molecule sets" 
[18-20]. 



Summarised result of the COSMO-RS calculations 

In aquo the hydrogen-bond model BCD023r061 was again the 
most preferred, not only among the empty (3-CD conformers but 
also as the BCD023r061/benzene parallel complex model. 
According to the quantum mechanical COSMO-RS calcula- 
tions, all empty (3-CD and all complex conformers of (3-CD/ 
benzene were thermodynamically allowed, i.e., the complexes 
should be formed. The highest relative energy of all conformers 
was less than 2.4 kcal/mol -1 in vacuo and less than 
6.4 kcal mol -1 in aquo. The benzene molecule adopted an 
oblique position in the inclusion complex. The o-surface of the 
empty (3-CD "looked" more hydrophilic from the outside at its 
02/03 side (red and blue) than from the 06 side (only red and 
yellow/green). From the inside, the 06 rim of (3-CD "looked" 
hydrophobic. 

Molecular-dynamics simulations 
The empty cyclodextrin at 300 K 

The "highly ordered hydrogen bonds" of the BCD0231061, 
BCD023r061, BCD023106r and BCD023r06r models cannot 
be observed at 300 K where the temperature is high enough to 
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Figure 6: a-Profiles of the COSMO-RS method for the four empty (3-CD models and benzene (BP/TZVP-DISP3 method). 
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Figure 7: o-Potentials of the COSMO-RS method for the four empty (3-CD models and benzene (BP/TZVP-DISP3 method). 
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make the P-CDs flexible and the energy barriers of hydrogen 
bonds are overcome. As was observed in the neutron diffrac- 
tion studies of P-CDs at 300 K [33] and 120 K [35] the tempera- 
ture determined the hydrogen-bond distribution. In the MD 
simulations started from the "highly ordered" 0 K AMI models 
the hydrogen bonds switched to more disordered arrangements. 
For comparison the crystal structure model of 300 K was added. 
The starting structures and the energetically best frame from the 



last 2000 ps of the 6000 ps trajectories are displayed in 
Figure 8, Figure 9 and Figure 10. Only the two energetically 
best models are shown: BCD023r061 and BCD023r06r. The 
Compass-istot energy of BCDcry was higher in the starting 
structure than in the best frame of the trajectory at 300 K, which 
accounts for the number of packing forces in the crystal that 
deformed the cyclodextrin monomer in vacuo. For the two other 
models BCD023r061 and BCD023r06r this was reversed; 




Figure 8: a: BCDcry (MD starting structure with Compass-E to t = 389.99 kcal mol -1 and 14 hydrogen bonds with averaged Dreiding-E H B = 

-3.04 kcal mol -1 => -42.56 kcal mol -1 ) and b: BCDcry fr1 662 (best frame from the last 2000 ps of the analysed MD-trajectory with Compass-E tot : 

308.01 kcal mol -1 and 11 hydrogen bonds with averaged Dreiding-E H B = -0.93 kcal mol -1 => -10.23 kcal mol -1 ). 




Figure 9: a: BCD023r06l (MD starting structure with Compass-E tot = 221.52 kcal mol -1 and 28 hydrogen bonds with averaged Dreiding-E H e = 
-4.71 kcal mol -1 => -131.88 kcal mol -1 ) and b: BCD023r06l fr4739 (best frame from the last 2000 ps of the analysed MD-trajectory with Compass- 
Etot = 292.47 kcal mol -1 and 13 hydrogen bonds with averaged Dreiding-E H B = -0.63 kcal mol -1 => -8.19 kcal mol -1 ). 
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Figure 10: a: BCD023r06r (MD starting structure with Compass-E to t = 217.25 kcal mol -1 and 28 hydrogen bonds with averaged Dreiding-Ene = 
-4.59 kcal mol -1 => -128.52 kcal mol -1 ) and b: BCD023r06r fr5418 (best frame from the last 2000 ps of the analysed MD -trajectory with Compass- 
Etot = 290.26 kcal mol -1 and 14 hydrogen bonds with averaged Dreiding-EnB = -0.65 kcal mol -1 => -9.1 kcal mol -1 ). 



their "highly ordered hydrogen-bond structures" in vacuo from 
the start were energetically lower than their structures from 
their best frames of the trajectories at 300 K, because they "lost 
their hydrogen-bond stabilisation" and became flexible instead. 
All starting structures had much higher hydrogen-bond ener- 
gies (analysis using Dreiding force field) than the best struc- 
tures from their trajectories. 

The five MD trajectories of the empty P-CD models had the 
following averaged energies (Compass force field): BCDcry 
311.15 kcal mol" 1 , BCD0231061 304.26 kcal mol" 1 , 
BCD023106r 304.02 kcal mol" 1 , BCD023r061 
295.14 kcal mol" 1 and BCD023r06r 291.22 kcal mol" 1 . 
According to this the BCD023r06r model was the lowest 
energy conformer, closely followed by the BCD023r061 
model, which was the best conformer in the AMI method. This 
means, that the hydrogen bonds at the larger 02/03 rim of 
empty P-CDs preferred the right-hand orientation, i.e., 
03-H 02-H in the same glucose unit and bifurcated towards 
04 and 03 of the next glucose unit on the right side. 

The p-CD/benzene complex at different tempera- 
tures from 300 K to 120 K 

The MD trajectories of the parallel and vertical BCD023r061 
model were practically the same since the benzene guest 
changed its position already in the first few pico seconds, and 
later a significant difference could not be observed any more. It 
was observed that the benzene guest left the cyclodextrin host at 
slightly different times, depending on the temperature, 
(Figure 11). At 120 K benzene was always inside the cyclodex- 



trin cavity. "Inside the cavity" was defined by a distance of up 
to 2 A between the centre of mass of benzene and the seven 04 
oxygen atoms of the P-CD that define a plane in the middle of 
the torus. At 273 K this distance extended up to 3 to 4 A 
starting at about 140 ps, and at about 320 ps it became longer 
than 5 A indicating that now benzene had left the cavity without 
return. At 280 K and 290 K these events appeared earlier, at 
140 ps/160 ps and 50 ps/125 ps, respectively. This was in good 
qualitatively agreement with the complex formation constants 
that were experimentally determined from UV spectra in [34] 
and showed that the complex was more stable at lower tempera- 
tures. At 300 K benzene tried to leave the cavity early at the 06 
side, but always returned. It spent most of the time at a greater 
distance of about 4 to 5 A to the 04-plane of P-CD and finally 
left the cavity on the 02/03 side as all the others did, but later 
at 2400 ps. Benzene had a high mobility in many geometrical 
positions in these trajectories, and the most characteristic ones 
from the trajectory at 300 K are shown in Figure 12. Figure 13 
displays how different the energies of the P-CD/benzene com- 
plex were during the MD trajectory at 290 K. 

Summarised result of the molecular dynamics simu- 
lations 

Molecular dynamics simulations at different temperatures 
(120 K to 300 K) displayed qualitatively correctly the dynamics 
of the P-CD/benzene complex. Benzene accepted many 
different positions inside the P-CD cavity before it finally left at 
the 02/03 side. The highly ordered hydrogen bonds at the 06 
or 02/03 rim were not existent at these temperatures and thus 
were a phenomenon of very low temperatures. 
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Figure 1 1 : Distance plots from the MD trajectories of 500 ps each {x axis) for (3-CD/benzene complex at different temperatures: a: 300 K, b: 290 K, 
c: 280 K, d: 273 K and e: 120 K. Distance [A] (y axis) between the centres of mass of benzene and the seven 04 oxygen atoms of [3-CD; horizontal 
red lines at distances of 2 A and 5 A, see text; vertical red lines indicate the time when the benzene left the (3-CD cavity without return. 
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Figure 12: Three structures from the MD trajectory at 300 K: 

a: benzene left the (3-CD at the 023 side; b: (3-CD/benzene complex; 

c: benzene tried to leave the (3-CD at the 06 side early, but always 

returned. Benzene finally left the cyclodextrin cavity also at the 023 

side. 



Since quantum-mechanical calculations lead to these highly 
ordered optimised structures, they represent the lowest energy 
conformers in vacuo at 0 K. Nevertheless, they will be the basis 
for the assumptions of the COSMO-RS method and its thermo- 
dynamic formulas and calculations/extrapolations based on the 
energetically lowest single-molecular states, plus conformers on 
top of that, up to 10 kcal mol -1 . These conformer ensembles 
will make it possible to calculate several thermodynamic prop- 
erties at room temperature and above, not only for the gas phase 
but also for the liquid and solid states of matter made from these 
molecules. 

Conclusion 

The P-CD/benzene AMI energies of complex formation after 
geometry optimisations for both positions of the guest, parallel 
and vertical were -2.75 kcal mol -1 and -2.46 kcal mol -1 , res- 
pectively. The concerted highly ordered hydrogen bonds at the 
02/03 and 06 rims of cyclodextrin had a strong influence on 
the structures. This was confirmed by the calculated AMI IR 
spectra, which showed that the P-CDs 0-H frequencies were 
mostly split up upon complex formation with benzene. The 
HOMO/LUMO gap of the empty cyclodextrin with about 12 eV 
was lowered to about 1 0 eV in the complex, in agreement with 
other calculations (PM3) from the literature. All four hydrogen- 
bond models studied here (BCD0231061, BCD023r061, 
BCD023106r and BCD023r06r) were energetically advanta- 
geous after geometry optimisations with BP86/TZVP-DISP3 in 
vacuo and in aquo, their relative energies were, all but one, less 
than 3.2 kcal mol -1 . The character of the COSMO-RS o-surface 
of P-CD was much more hydrophobic on its 06 rim than on its 
023 side when all hydrogen bonds were arranged in a concerted 
mode. COSMO-RS BP/TZVP-DISP3 calculations in vacuo and 
in aquo (dielectric field approximation) preferred an oblique 
position (over parallel or vertical) of benzene inside the P-CD 
cavity and suggested energies of complex formation up to 
20.2 kcal mol -1 in vacuo and of about 16 kcal/mol -1 in aquo. 
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The static picture of these energetically lowest structures of the 
P-CD/benzene complex at 0 K was extended by the molecular 
dynamics simulations at different temperatures (120 K to 
300 K). They displayed qualitatively correctly the dynamics of 
the P-CD/benzene complex, which was more stable at lower 
temperatures, especially at 120 K. Benzene accepted many 
different positions inside the P-CD cavity before it finally left at 
the 02/03 side (273 to 300 K). The highly ordered hydrogen 
bonds at the 06 or 02/03 rim were a subtle phenomenon at 
very low temperatures. Still, we assume that hydrogen bonds of 
cyclodextrins nevertheless play a crucial role in all their special 
properties, i.e., their many inclusion complexes and remarkable 
influence on reactions, catalysis and supramolecular structures. 
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